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The total momentum of quantum fluids 
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The probability distribution of the total momentum P is studied in A-particle interacting homoge¬ 
neous quantum systems at positive temperatures. Using Galilean invariance we prove that in one 
dimension the asymptotic distribution of P/ y/N is normal at all temperatures and densities, and in 
two dimensions the tail distribution of P/ y/N is normal. We introduce the notion of the density ma¬ 
trix reduced to the center of mass, and show that its eigenvalues are N times the probabilities of the 
different eigenvalues of P. A series of results is presented for the limit of sequences of positive defi¬ 
nite atomic probability measures, relevant for the probability distribution of both the single-particle 
and the total momentum. The P = 0 ensemble is shown to be equivalent to the canonical ensemble. 
Through some conjectures we associate the properties of the asymptotic distribution of the total 
momentum with the characteristics of fluid, solid, and superfluid phases. Our main suggestion is 
that in interacting quantum systems above one dimension, in infinite space, the total momentum is 
finite with a nonzero probability at all temperatures and densities. In solids this probability is 1, 
and in a crystal it is distributed on a lattice. Since it is less than 1 in two dimensions, we conclude 
that a 2D system is always in a fluid phase; so if the hexatic phase existed classically, it would be 
destroyed by quantum fluctuations. For a superfluid we conjecture that the total momentum is zero 
with a nonzero probability and otherwise its distribution is continuous. We define a macroscopic 
wave function based on the density matrix reduced to the center of mass. We discuss how dissipation 
can give rise to a critical velocity, predict the temperature dependence of the latter, and prove that 
Landau’s criterion cannot explain superfluidity and its breakdown in a dissipative flow. We also 
comment on the relation between superfluidity and Bose-Einstein condensation. 
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I. INTRODUCTION 

Proving phase transitions in continuous space is notori¬ 
ously difficult. To the author’s knowledge there are only 
three rigorous results for classical systems, by Ruelle [1] 
about segregation in the Widom-Rowlinson model, by 
Lebowitz, Mazel and Presutti [2] about the vapor-liquid 
transition, and by Bowen et al. [3] about fluid-solid tran¬ 
sition, and none for interacting homogeneous quantum 
systems. Here we suggest the use of the total momentum 
for the characterization of the different thermodynamic 
phases of quantum systems in continuous space, at posi¬ 
tive temperatures. A short account of this work is given 
elsewhere [4]. 

At a first sight the total momentum is not a very excit¬ 
ing quantity. In finite volumes it is conserved only if the 
boundary condition is periodic. Moreover, because it is 
the sum of N single-particle momenta, one would guess 
it to be typically of order y/lv. Yet a second thought may 
modify this picture. Thermodynamic phases can unam¬ 
biguously be characterized only in the thermodynamic 
limit, so if we are interested in homogeneous systems, 
the use of periodic boundary conditions is not really re¬ 
strictive. More important, the ‘typical’ y/N behavior 
does not mean exclusivity. Indeed, in interacting quan¬ 
tum systems there are excitations in which a macroscopic 
number of particles is involved. If these ‘collective modes’ 
carry a momentum, as in the case of density waves and 
lattice vibrations, this is, or is proportional to, the total 
momentum; so manifestly the total momentum can re¬ 
main finite in the thermodynamic limit. The center of 


our interest is the probability distribution of the total 
momentum in thermal equilibrium. 

In Section H we write down a general N-particle 
Hamiltonian, the corresponding density matrix and par¬ 
tition function in d-dimensional cubes with periodic 
boundary conditions. We then use Galilean invariance 
to sum separately over eigenstates belonging to the ir¬ 
reducible wave vectors, introduced in Ref. [5], and over 
their Galilean boosts. This leads to some interesting re¬ 
sults. The set of irreducible wave vectors is markedly 
different in one, two, and higher dimensions. In one di¬ 
mension we can derive a central limit theorem showing 
that the limit of Q/y/lV, where HQ are the eigenvalues 
of the total momentum operator, is normally distributed 
at all temperatures and densities (Theorem H.I). In two 
dimensions only the tail distribution is normal (Theorem 
H.2), and for higher dimensions we can prove only that 
IQI/iV = 0{1/L) almost surely (Proposition H.3). 

In two theorems of Section HI we investigate the rela¬ 
tion between the total momentum and the center of mass. 
In Theorem HI.I first we show that in any eigenstate of 
the total momentum operator the center of mass can be 
separated and it propagates as a free particle carrying 
the total momentum. This result, applied already in our 
earlier paper [5], is intuitively obvious but its derivation 
needs some care. We then use this separability to define 
Pc.m., the density matrix reduced to the center of mass, 
and prove that the eigenvalues of Pc.m. are N times the 
probabilities of the different eigenvalues of the total mo¬ 
mentum. The integral kernel (xjpc.m. |y) is shown to be 
real nonnegative; we interpret it as the autocorrelation 
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function of the center of mass. Theorem III.2 extends the 
definition of Pc.m. to cases when a restricted density ma¬ 
trix is reduced to the center of mass; e.g. this can be the 
projection of the density matrix to the Q = 0 subspace. 

Section IV is a brief summary of some known results for 
the single-particle momentum distribution. We include it 
to stress the structural analogy between pc.m. and pi, the 
one-particle reduced density matrix. The expected frac¬ 
tion of particles in the one-particle state of momentum 
fik is the analog of the probability that the total momen¬ 
tum is /ik, and the integral kernel (x|pi|y) is also real 
nonnegative; moreover, (x|pc.m.|x) = (x|pi|x) = p, the 
global density. The fundamental difference between Pc.m. 
and Pi is due to the fact that the single-particle motion is 
not separable [except for the noninteracting Boltzmann 
gas]. 

Section V is about the infinite volume limit of the prob¬ 
ability distributions introduced for the total momentum 
in Section II and for the single-particle momentum in 
Section IV. In both cases we are given a sequence ipL 
of positive dehnite atomic probability measures concen¬ 
trated on a lattice whose spacing tends to zero as l/L. 
The distribution limit of LpL is a not necessarily normal¬ 
ized positive measure in the continuous dual space. We 
can compactify the latter by adding a point at infinity, 
and attribute the missing weight to that point. The sec¬ 
tion contains a series of propositions and theorems relat¬ 
ing the weight at infinity to the properties of the limit 
of (pz, and of its Fourier transform. We also give exam¬ 
ples that we think to be characteristic to crystals, fluids, 
Bose-condensates, and superfluids. 

In Section VI we show that from the point of view 
of thermodynamic quantities the subspace of zero total 
momentum is representative to the whole system: the 
asymptotic free energy density can be computed in the 
Q = 0 restricted ensemble. With this restriction we lose 
information about the probability distribution of Q, with 
a notable exception: Theorem III.2 implies that the prob¬ 
ability of Q = 0 in the canonical ensemble can directly 
be obtained from the reduction of the density matrix to 
the center of mass within the Q = 0 subspace. 

In Section VII we make some conjectures about the 
relation between the probability distribution of the total 
momentum and the type of the thermodynamic phase, 
partly anticipated in the examples of Section V. The to¬ 
tal momentum is the property of a state of the whole 
system. So when we say that in an infinite fluid (gas 
or liquid) the total momentum can be both finite and 
infinite, we mean that the thermal equilibrium state is 
a convex combination of states with different total mo¬ 
menta, both finite and infinite. We conjecture that in an 
infinite solid the random few-particle motion freezes out, 
the total momentum is finite with probability one and has 
a purely atomic probability distribution; in crystals the 
probability distribution is concentrated on a lattice. [In a 
supersolid one would have the convex combination of an 
atomic and a continuous probability distribution.] This 
corresponds to what we found in one dimension at zero 


temperature [5]. Because we proved earlier that in two 
dimensions at all temperatures and densities the total 
momentum diverges as with a nonvanishing proba¬ 
bility, cf. Theorem II.2, we conclude that no solid phase 
can exist in 2D; so quantum fluctuations would destroy 
the hexatic phase in case it existed classically. 

The most audacious among the conjectures is the one 
about the superfluid. Superfluidity is well understood ex¬ 
perimentally, and there exists a fully developed hydrody¬ 
namic theory initiated by London [7], Tisza [8] and Lan¬ 
dau [6] which, assisted by the weaponry of many-body 
physics, provides a satisfactory explanation to practically 
all experimental findings. Moreover, this theory is acces¬ 
sible in monographs written by top experts of the sub¬ 
ject [7, 9-12]. However, the first-principle characteriza¬ 
tion of superfluids is a slippery problem [13]. A pioneer¬ 
ing work about the A-transition in liquid helium, based 
on first principles, is due to Feynman [14]. Through his 
path integral representation of the partition function, by 
making a series of approximations he obtained a third 
order phase transition with a continuously varying spe¬ 
cific heat, instead of a discontinuous one. His method 
has remained the most promising route to a proof, but 
until today nobody was able to accomplish the task. 

Our conjecture says that below the superfluid tran¬ 
sition temperature the total momentum of the infinite 
system is zero with a nonzero probability, and the state 
of zero momentum can be identified with the superfluid 
component. The normal component is the state in which 
the total momentum is nonzero, and the thermal equi¬ 
librium state of the infinite system is the convex combi¬ 
nation of the two components taken with the respective 
probabilities. A measurement can project the system into 
one or the other component. This definition makes the 
superfluid transition reminiscent to Bose-Einstein con¬ 
densation (BEC), only the ’condensation’ takes place into 
the state of zero total [and not single-particle] momen¬ 
tum. The discussion of the consequences of this conjec¬ 
ture occupies the major part of Section VH. Our defini¬ 
tion is motivated partly by aesthetic reasons. First, there 
is the appealing analogy between pc.m. and pi. Second, 
as the temperature goes to zero, the Gibbs state tends 
to the ground state which is both 100% superfluid and 
of zero total momentum. Third, it follows that in the 
superfluid state thus defined the system really moves ‘as 
a whole’, as Landau phrased it [6]. At rest or in a hy¬ 
pothetic frictionless flow obtained by Galilean boost the 
superfluid is in thermal equilibrium; in a realistic dis¬ 
sipative flow the superfluid is a non-equilibrium state. 
The effect of dissipation can be understood as the break¬ 
down of thermal equilibrium in such a way that in finite 
volumes the density matrix becomes the convex combi¬ 
nation of two states of equal temperatures but generated 
by different Hamiltonians. The existence of a critical ve¬ 
locity together with its temperature dependence ensues 
from this picture. A different view of the non-equilibrium 
nature of the superfluid flow at zero temperature is pre¬ 
sented in Ref. [15]. 



3 


At the end of Section VII we comment on the rela¬ 
tion between superfluidity and Bose-Einstein condensa¬ 
tion. In two dimensions there is no BEC at positive tem¬ 
peratures [16-18], but a superfluid transition takes place 
in trapped 2D Bose gases [19] and in helium monolay¬ 
ers [20]. Bulk liquid helium undergoes both BEC and 
superfluid phase transition, and there is today a consen¬ 
sus that Fritz London’s original idea [7] was correct, and 
the superfluid phase transition and BEC occur simultane¬ 
ously. Yet even in 3D one cannot equate the superfluid 
with the Bose condensate: as the temperature goes to 
zero, the whole system becomes superfluid while the Bose 
condensate contains less than 10% of the particles [21- 
26]. The condensate wave function, obtained from pi, is 
generally accepted to be the macroscopic wave function 
associated with the superfluid. However, a macroscopic 
wave function derived from pc.m. can play the same role 
even better: because the superfluid fraction, as we define 
it, tends correctly to 1 with the temperature going to 
zero, one can avoid the embarrassing question, how less 
than 10% can represent 100%. Meanwhile, we do not 
challenge the consensus. In three dimensions BEC and 
the superfluid transition may well be simultaneous. Also, 
there is nothing contradictory in supposing that the total 
momentum and a macroscopic number of single-particle 
momenta vanish simultaneously. The advantage of our 
proposal is that it points to a common cause of superflu¬ 
idity in two and three dimensions. 

In Section VIII we discuss Landau’s criterion for the 
breakdown of superfluidity. Doubts about its validity 
were raised already earlier, see e.g. in Refs. [13, 27]. Here 
we analyze Landau’s original publication [6] and find that 
his argument does not apply to a dissipative flow. The 
paper ends with a Summary. 

II. DISTRIBUTION OF THE TOTAL MOMEN¬ 
TUM 

The setting is the same as in Ref. [5]. 

1 ^ 

+^A(xi,...,XAr) (1) 

is the energy operator of N interacting particles in a el- 
dimensional cube A of side length L, defined with peri¬ 
odic boundary conditions; Ua is the periodized potential 
energy. Stability is supposed, and hard-core interactions 
are allowed, with the suitable modification of the domain 
of H. H is either unrestricted [Boltzmann statistics] or 
restricted to the symmetric or antisymmetric subspace. 
The common eigenstates of H and the total momentum 
operator 


The energies are ordered, < Uq_„+i (n > 0). The 
ground state and its energy are V’o.o and Uq.o, respec¬ 
tively. With the notation 

IIq,„ = |V’Q,n)(V’Q,n| (4) 

the unnormalized density matrix is 

OO 

E (5) 

QgA* n—0 


Introducing the set of irreducible wave vectors 

^kred = {q G A* : —ttN/L < Qi < ttN/L all i}, (6) 

any Q G A* can uniquely be written as q -I- Vk, where 
q G At and k G A*. Then, 

OO 

E E (7) 

qeAfrred k6A* "=0 
By Galilean invariance [5], 

Ec^+Nk,n = £^q,n + /2m) [Ak^ -h 2k • q]. (8) 

Substituting this into Eq. (7), 

OO 

^ ^ g-/3Sc,,n 

qeVrred ”=« 

X ^ (9) 

keA* 

In what follows, we fix the number density p; then 
N = pL’^. With A /3 = A/27r/3fi2/m, the partition function 
reads 


Z= Tre"^^= ^ ^g-/5£;c,,n ^ g 

qeVrred ”=0 kSA' 


Let 


-^irred — E E^ 

q6A*rred "=0 


-/SiJq,. 


-if [AfkU2k-q]^ 

( 10 ) 

( 11 ) 


the k = 0 term of Z. Then ZjZA-ced is the average in 
At of the sum over k. 


Zirred 




( 12 ) 


irred 


p = EPi = 

J 3 


( 2 ) 


First we bound this quantity. 

Proposition II.1. 


are denoted by ipQ.n', the respective eigenvalues are Uq^n, 
and hQ, where 


(3) 


max 





< 


-^irred 


< 


r^l-d/2 n 

L J 


Q G A* = 


(13) 
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Proof. The sum over k in Eq. (10) can be factorized. 
Replacing the ith component of k by 271 j /L, 


H e 4 
keA’ 


^[iVk=+2k.q] _ 


n 


E 


e ^ cosh 


n 


1 + 2E 

i=i 


e ^ cosh 


L 


By taking the limit L —>■ oo, a first estimate about the 
asymptotic behavior of E^ can be obtained. 

Proposition II. 2. For any /3 and p, 

• if d = 1, 


lim {Tip) = 0; 

L—¥(yo 


(14) 


• if d = 2, 


For q G and j > 0, 


1 / k / 7rA|lVj 

1 < cosh — < exp ■ 


(15) 


L - L2 ■ 

Equation (13) follows from the upper and lower bounds 


E^ ^ 

1=1 


\j,N 

—f4-Aj-i) 


< 1 


^OO nXiN 
- 7 ^—a 

e 


dx 


and 


e r-4 i > jnax "I 0, 

1=1 


e L 


2 

dx-1 


(16) 


(17) 


with 


e L 


2 

“T-^ 


pi-dn 

dx = — -—. □ 

2A/?vr 


Let 






,-/3Bc,,„„-^[lVk2+2k.q] 


n —0 


n=0 


the probability that in thermal equilibrium the system 
is in the state with fiQ = 

ft(q + iVk). By definition, X^qeA* = 1) meaning that 
in finite volumes the total momentum is finite with prob¬ 
ability 1. As we shall see, this is not the case in infinite 
volume. Let us introduce also the probability distribu¬ 
tion function 


Tl (k) = E '"Q 

IIQII<« 


(19) 


where ||Q|| = maxi<j<d jQ^j. In words, r^(«;) is the 
probability that the modulus of each component of the 
total momentum is less than or equal to Hk. Then 


■^irred 


= E ^Q^^Li^N/L)., 


( 20 ) 


cf. Eq. (6), and Eq. (13) can be put in the form 


< X! < min <1 1, 


L 


-f 3 


qGA* 


Ll~d/2 


1 d 


^fiVp 


- 1 


( 21 ) 


lim i'^pL) > 




L —¥00 


{\ + d\p^y 


if d > 3, 


lim {TrpL’^~y > l/S”^. 

L—>oo ' ' 


( 22 ) 


(23) 


(24) 


Proof. For d = 1 the upper bound, for d > 2 the lower 
bound in Eq. (21) provides the result. □ 

Because Q is the sum of N single-particle wave vec¬ 
tors, normal fluctuations of the total momentum about 
the zero expectation value would mean that the absolute 
value of each component of Q is typically of order VN. 
Let J be any positive integer and consider the probability 


'' ^ |Q^|> (2J-0.iV 


l^Q. 


From Eqs. (14) and (18), 


(25) 


( 18 ) F2J-1 = ^E E 




n=0 qGA* 


E 


e ^ cosh — - 


d-l 

n 

2=2 


> e ^ cosh 




J^-OO 


(26) 


Theorem II. 1. In one dimension the asymptotic dis¬ 
tribution of Q /^/N is normal at all temperatures and 
densities. Namely, 

lim Prob I < xl = dy. (27) 

\ p^/N ~ j tt 7o 

Proof. From Eq. (26) we subtract the same equation 
written for J -|- I: 

2 “ „ _ X^qj 

F2J-1-F2J+1 =—^ E ® ^ ^ cosh^^. 

n =0 96A*„^_j 

(28) 

Because |( 7 | < vrp if g G summing over J from 1 to 

M and using the bounds 

I < cosh [yf^qJ/L) < cosh (7rA|pAf/L) (29) 
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we find 

< F, - F^m+1 

j=i 

< cosh {TrXjpMlL) ^ ( 30 ) 

J=i 

Now 

Fi - F 2 M +1 = [(2M + l)7rp] - (irp), (31) 

where r^(7rp) —>• 0, c.f. Eq. (22). Setting M = x'/N, 
cosh ^TrA^pM/L^ —>■ 1. Furthermore, 

1 xsTN „x 

lim ^ V dy. (32) 

With the bounds (21), 

= 2\pp [ dy 

Jo 

(33) 

which is equivalent to Eq. (27). □ 

The above result obviously implies that for any se¬ 
quence ol —>■ 00 , 


lim r 
N—¥oo 


N 


{2xVn + l)7rp 


Proof. With J = 2 the first sum over j in Eq. (26) goes 
to zero while the remaining sums go to 1. So Ea —>■ 0 
which gives the result. □ 

Above two dimensions we do not obtain the expected 
y/N instead of because for d > 3 any Q of or¬ 
der VN is deeply inside and one should conceive 

a proof which takes the interaction into account. An 
idea for such a proof is as follows. The result for one 
and two dimensions implies that the thermal average 
(p2) = 0{N). Because (p^) is finite, (Pi-Pj) = 0(1/A^). 
The argument can work in the other direction. If one 
could prove that (pi • pj) = 0(l/iV), one could conclude 
that (P^) = 0{N). If the scattering length of the inter¬ 
action is finite, the momentum of each particle is corre¬ 
lated only with that of other particles within scattering 
distance. But for fixed x^, Xj can be anywhere in A, 
therefore (pi • p^) = 0{L~‘^) = 0(1/N). Note that ir¬ 
respective of an eventual improvement of the result for 
d > 3, in all dimensions the velocity P/Nm of the cen¬ 
ter of mass tends to zero with probability 1 as L goes to 
infinity. 

III. TOTAL MOMENTUM AND CENTER OF 
MASS 

Let X denote the arithmetic mean (center of mass) of 
the vectors xi,..., x^r. It is easy to check that for any 
differentiable function /(xi,..., x^v). 


\im (p/Aol) =1, lim T^ (^p^N/ol) = 0. 

(34) 

Theorem II. 2. In two dimensions the asymptotic tail 
distribution of |Qi|/-\/]V is normal at all temperatures 
and densities. More precisely, for any J > 1, 

2 /(A^v^!) Prob{JM>(2J-l 

j=J C V P 

00 

j=J 


(x • P — P • x)/ = idhf, (38) 

that is, the total momentum is canonically conjugate to 
the center of mass. A more interesting relation between 
the two is as follows. 

Theorem III.l. {A^^g,are the eigen¬ 
values and corresponding eigenfunctions of a reduced 
Mensity matrix, where the reduction is onto the center 
^ of mass. In Boltzmann and Bose statistics, vq is a pos¬ 
itive definite function of Q on A*, i.e., the n x n ma- 
(35)trix is positive semidefinite for any n and any 

Qi,..., Q„ e A*. In particular, 


where 


! irred 

(36) 

Proof. We apply Eq. (26) with the lower and upper 
bounds obtained from Eq. (15). / (A^^p) > 0 follows 
from Eq. (23). Xp,/p is the only relevant dimensionless 
quantity, whence the dependence of / on it. For J large 
enough the upper bound in Eq. (35) is less than 1. Thus, 
the limit of the probability is between 0 and 1 and tends 
to zero as with J going to infinity. □ 

Equations (21) and (26) have only a much weaker im¬ 
plication for d> 3. 

Proposition II.3. If d > 3 then for any /3 and p. 


lim rf (37rA/L) = 1. 

L—^oo 


> vq, any Q G A*. (39) 

Proof, (i) Separability of the center of mass. Any eigen¬ 
state of the total momentum, hence, also V'Q.n can be 
written as 


V'Q.n(xi, . . . ,X]v) = 


(40) 


where 


N 


= > X)v), X' = Xj - 


Xl 


i=i 


(41) 


[so = 0]. Indeed, for any a G 

e ^ ^^/iQ^„(xi,... ,XAr) = e“*‘^'‘^'i/'Q.«(xi,... ,xjv), 

(42) 


(37) 
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and also 

e " ... ,X7v) = V'Q,«(xi - a,... ,xjv - a). 

(43) 

Comparing the two, 

V’Q,n (xi,..., Xiv) = (xi - a,..., XAT - a). (44) 

From Eq. (44) with a = xi and 

N 

xi=x —7V“^^x' (45) 

i=2 

we find 

V>Q,„(xi,. ..,^n) = ^^qAO, X'). 

(46) 

Introduce 

X = (xi,... ,XAr), Y = (yi,...,yAr) (47) 

with 


Fourier transformation of Eq. (52) in x and y shows that 
Pc.m. is diagonal in momentum representation, 

(QIPc.m.lQO =(54) 

so its eigenvalues are Ni^q, its spectral resolution reads 
Pc.m. =iV ^ J^q|Q)(Q| (55) 

QeA* 

with (x|Q) = 

(ii) Positive definiteness of vq. We prove that 

(x|pc.m.|y) > 0. (56) 

Pc.m. can also be obtained from the integral kernel of 
, without passing through the spectral resolution of 
the latter. With the same X and Y as above, in the 
Boltzmann case we obtain 

(x|pe.m.|y) = f / dx' (X|e-^^|Y) ; (57) 


X, = x + x' - ^ ^x', y, = y + x' - ^ ^x', (48) 

i—2 i—2 

where x, y, X 2 ,..., x(y are arbitrary and x( = 0. Then 
Eq. (41) holds true, so x is the center of mass of {xj}, 
and one can check that y is the center of mass of {yj}. 
Equation (46) is valid and, because yj — Xj = y — x 
independently of j, one can apply Eq. (44) with a = x—y, 
yielding 

V^q.„(Y) = 

= e*‘4'ye-(*/^)‘^'Xf=2>^;q/,Q_„(0,X'). (49) 

As a result, 

V-Q.„(X)V>^_„(Y) = e*Q-(x-y) ^(0, x')|' • (50) 

Because of the periodic boundary condition, integration 
over X' in A^v- 

gives the same result as integration over 
X 2 ,..., xvr. Therefore, using again Eq. (44), 

J |V>Q.„(0,X')|" dX'= J |V>Q.„(X)|" dx 2 ---dxvr = T-'' 

(51) 

independently of xi. Now we can define Pc.m., the den¬ 
sity matrix reduced to the center of mass [28], through a 
partial trace over X': 

(Xjpe.m.ly) = iVZ-l^e-^^«-/'V-Q.n(X)V^^_„(Y)dX' 

Q,n 

= d E (52) 

qga* 

It is seen that 

(xjpc.m.ly) < (x|pc.m.|x) = P E = P' (53) 

qga* 


for Bose statistics 

(x|p.,„,|y) = 5j^ E / dX'(GX|e--»|Y) 

GgSn 

(58) 

where S'at is the symmetric group of N elements, and 

GX = (xg(i), •. • ,xg(jv)) (59) 

[here H is unrestricted; restriction to the symmetric sub¬ 
space is explicitly taken into account]. Equation (56) 
follows from 

(Xle-'^^IY) > 0 (60) 


which can be traced back through the Trotter formula to 


(x|e-“'P.^|y)=L-^E 


.-“|k|^pik-(x-y) 


e "‘o 


e-^^l"l"-^"4x-y) > 0, (61) 


obtained via the Poisson summation formula [a = 
a'h^/2m]. Therefore, 

Nvq = (QjPc.m.lQ) = ^^^(x|Pc.m.|y)e*‘^'^^"’"Mxdy 

= /(0|pc.m.|z)e*‘^ '^dz < (62) 

JA 

Being the Fourier transform of an integrable positive 
function, vq is a positive definite function (a function 
of positive type [29]). □ 

Because > 0, by Eqs. (55) and (52) pc.m. is a 
positive operator, and its kernel is also a positive defi¬ 
nite function: the n x n matrix [(xi|pc.m.|xj)]ij is posi¬ 
tive semidefinite for any n and Xi,..., x„ G A. Clearly, 



7 


(x|pc.m.|y) = (y|Pc.m.|x) Or (0|pc.m.|z) = (0|pc.m.|- z) • 
One can interpret p“^(x|/9c.m.|y) as the autocorrelation 
function of the center of mass. 

Note that Eq. (60) remains valid in the subspace of 
zero total momentum. We prove a more general result. 
Introduce 


OO 

nQ = ^nQ,„. (63) 

n—0 

By definition (4), X^qeA* = In, the identity operator 
in 2.2 (A^). 

Theorem III. 2. Let L* be any subgroup of A*, and let 
4> = {'/’q} be a positive definite function on L*. Then 




-0H 


qeL* 


Y ) > 0. 


(64) 


Moreover, if we define by its integral kernel as 


mptm.rn = Nz-^ J (: 


'/’Qnq' 


-/Sff 




where X and Y are given by Eqs. (47), (48), then 

(x|Pc.m.|y) = P > 0 . 

qeL* 


Yj dX' 
(65) 


( 66 ) 


In particular, with (()q = dq o we obtain 

(X|noe-^"| Y) > 0 (67) 

and Ilopc.m. = A^i^o|0)(0| or 

(x|noPc.m.|y) = pt'o- (68) 

Proof. Let X — a = (xi — a,..., xjv — a). Then 

(X|nqe-^^|Y) = ^ (X) V’q,„ (Y) 

n 

= E <^Q',QE®”''''‘"''’‘V'Q'.«(X)V’Qq„(Y) 
q'eA* n 

= L-^ [ dae*Q-^e-^^Q'."V’q'.«(X-a)V>a!y„(Y) 

q',n 

= L-^ [ dae*'5 ‘"(X-a|e-^^|Y) . (69) 

JA 

Multiplying by ^q and summing over Q in L* we find 


Reduction of Ilqe to the center of mass provides di¬ 
rectly J^q|Q)(Q|, from which Eq. (66) follows. □ 

IV. DISTRIBUTION OF THE SINGLE¬ 
PARTICLE MOMENTUM 

For comparison, here we collect some known results. 
The occupation-number operator Ak for the one-particle 
state |k) [k e A*] is a sum of projectors, 

N 

= J = Ij-i ® |k)(k| 0 In-j . (70) 

i=i 

Ak is not conserved; we have access to its thermal ex¬ 
pectation value (Ak) through the one-particle reduced 
density matrix, of which it is the eigenvalue [21]. Indeed, 
for bosons 

(Ak) = Z-^Tt [e-^^Nk] 

= E (ki, ■ • ■, kjvj Akue^^^Akujki,..., kjv) 

ki,...,kjv 

= ^ E (k,k2,...,kAr|e"'^^|k,k2,...,k7v) 

k2,...,kAr 

= (k|di|k), (71) 

where 

pi = AZ-iTr2....,jve-'^^ (72) 

is the one-particle reduced density matrix. However, pi 
is diagonal in momentum representation: because H con¬ 
serves the total momentum, 

(k|pi|k') = ^ E (k,k2,...,kAr|e"'^^|k',k2,...,kAr) 
k2,...,kA2 

= dk.k'(k|pi|k). (73) 

Therefore, 


di = E (^k)|k)(k|, (74) 

keA* 

and its integral kernel is 

0 < (xjpijy) = L-^ ^ (Ak)e*‘"'(’'~y) < (xjpijx) = p. 
keA* 

(75) 

The kernel is nonnegative because it can also be obtained 
by integration from the integral kernel of . Now 

0 < rik = (Ak)/A <1, ^ uk = 1, (76) 

keA* 


E > 0 

qeL* 

(in distribution sense) which together with Eq. (60) gives 
Eq. (64). Above we used Eq. (44). It is seen that 
(X|nqe“^^|Y) is also a positive definite function of Q. 


therefore nk can be interpreted as the probability of the 
occurrence of a particle in the one-particle state |k). We 
can also define 

ri(«)= E (77) 

kG A*: ||k|| <K 



the probability distribution function for the single¬ 
particle momentum. Altogether, the conclusion is the 
same as for pc.m.- 

Proposition IV. 1. {Vnk, are the eigen¬ 

values and corresponding eigenfunctions of the one- 
particle reduced density matrix pi. The integral kernel 
(0|pi|z) is a nonnegative positive definite function of z 
on Moreover, rzk is a positive definite function of k 
on A*; in particular, no > nk- 

Thus, Pi has the same structural properties as Pc.m., 
and nk is the analogue of nk, k going over the same set 
A*. However, their k-dependence or the ^-dependence 
of and is expected to be quite different, owing to 
the fact that the dependence of i/'Q,n(xi,..., x^r) on x is 
separable while its dependence on Xj is not. In what fol¬ 
lows, we investigate the possible alternatives, with special 
attention to the limit of infinite space. 

V. LIMIT OF SEQUENCES OF POSITIVE 
DEFINITE ATOMIC PROBABILITY MEA¬ 
SURES 

In both cases, that of the total and of the single¬ 
particle momentum, we are given a sequence of cubes 
A = [—L/2,L/2]‘^ whose side length L tends to infinity, 
the associated dual lattices A* = (27r/L)Z‘^, and for each 
L an even probability measure pL on A*: 

(/?L(k) = (/?L(-k) > 0, ^(/?z,(k) = l. (78) 

keA* 

Thus, the Fourier transform 

A(x) = ^ (79) 

keA* 

is a real even continuous positive definite A-periodic func¬ 
tion on fL{^) < /l( 0) = 1- Moreover, we know that 
fL{^) > 0. This implies that 

PL{k)=L-^ [ (80) 

JA 

is positive definite on A*. In the present applications 
(/?L(k) is either j^k or nk and, accordingly, /l(x) is either 
P“^(0|pc.m.|x) or p“1(0|/9i|x). 

For K > 0 let 

rL(K) = 

kG A* ,||k|| </* 

and 

r(n) = limsup ri(K). (82) 

L—¥00 

The corresponding quantity for the total and single¬ 
particle momenta are T'^ ™' and F^, respectively. 

Proposition V.l. r(K) is monotone increasing. 

Proof. Let k < k', and let Li be a subsequence such that 

lim rLi(n) = r(K). 

i—^oo 


Because is monotone increasing for any L, 
rL,(n)<Fi,(n') 

for each i. By the Bolzano-Weierstrass theorem, there is 
a subsequence Li^ of Li on which convergence takes place 
at k' , and necessarily 

F(k') = limsupFi(K') > lim F^^ («;') > F(a). □ 

L->oo ^ 

To discuss limits when L goes to infinity, it is useful to 
extend pL from A* to (R'^)* by using Dirac deltas: 

PL= ^L(k) = V3L(q)5(k-q). 

qSA* qGA* 

(83) 

Then 

/l(x) = y (^z,(k)e"**" ’'dk (84) 

and, if k is not an integer multiple of 27r/L then 

rL(K) = [ (^L(k)dk. (85) 

d||k||<« 

li Pl converges to p in distribution sense then F/, is 
convergent and 

r('«) = [ (^(k) dk. (86) 

J\M<^ 

This holds because the characteristic function of the cube 
{||k|| < k} can be approximated by functions of rapid 
decrease. However, p may not be a probability measure, 
because 


F<oo = lim F(k) = [ p{k)dk 

K—^OO J 

(87) 

may be less than I. We can remedy this by a one-point 
compactification of (R'^)*, 

(R^)* = (R^)* U {oo}. 

(88) 

and defining 


Foo = I ~ F<oo. 

(89) 

Let 


Ga(x) = (7rA")-‘^/2g-xVA^^ 

(90) 


When A goes to zero, G\ tends to (5o in distribution sense. 

Proposition V.2. If p^ tends to p in distribution sense 
then 

/(x) := lim lim {fL*Gx){x) = [ dk(^(k)e“**"’", (91) 

A—>0 L—>-oo J 

implying the existence of the limit, and 


r<oo = /(o). 


(92) 
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Proof. 

(/l*Ga)(x) = J dk^L{k)e-^'^-^ I dyP^-yGx{y) 

= J dk(^z,(k)e-**" ’'e-^'‘""/4. (93) 

Because in the integrand f>L{k) is multiplied with a func¬ 
tion of rapid decrease, one has 

lim (/l * Ga)(x) = [ (94) 

L^oo J 

Setting A = 0 and x = 0, comparison with Eq. (87) yields 
the result. □ 

We define 

/(x) = lim /l(x) (95) 

L—^oo 

provided that the limit exists. If this is the case, the 
question naturally arises whether f = f. If Too > 0, the 
answer is not, because 

/(0)-/(0) = i-r<oo = roo. (96) 

However, 

/(x) = /(x) if X ^ 0 (97) 

may hold true. Actually, this can easily be seen in one 
dimension, because for x ^ 0, coskx can be piecewise 
approximated by functions of rapid decrease. In the ex¬ 
amples given below we shall assume Eq. (97) also for 
d > 1. This implies that 

r<oo = lim/(x), (98) 

x->0 

so that if r<oo < I then /(x) is discontinuous at zero 
with a jump Too. 

r<oo and Too can be interpreted as the probabilities 
that in the limit of infinite space k is finite or infinite, 
respectively. Here is a quick result about the positivity 
of Too. 

Proposition V.3. Suppose that there is a sequence kl 
such that limi_>oo kl = co and 

limsupri(KL) < 1. 

L—^OO 

Then Too > 0. 

Proof. For any k < oo, rL(K) < Tl{kl) if L is large 
enough, and thus 

r(K) < limsupri(Ki). 

L—¥00 

Because the right member is independent of k, we can 
take the limit k —>■ oo and obtain r<oo < I. □ 

Corollary. In one dimension, applying the proposition 
to limsup^_,,gQ with kl = P\fpL/aL where ol —>■ 


oo, gl = o{'/L), from Eq. (34) we obtain that for all 
/3 and p, = 1; in two dimensions, from Eq. (35) 

with large enough J it follows that > 0. Later we 
shall associate these findings with the absence of phase 
transitions in ID and the absence of fluid-solid phase 
transition in 2D. 

Examples, (i) Let bi,...,bci be linearly independent 
vectors and 

L* = {27rnibi H-h 2T:ndhd}m,-,n,i&- (99) 

Suppose that the extension of to (R"^)* tends in dis¬ 
tribution sense to an atomic probability measure concen¬ 
trated on L*: 

f'L = ^ ^ (100) 

k^A* 

Then the infinite space limit of the autocorrelation func¬ 
tion of the center of mass, 

/'=■“■ (x) = lim p"^(0|pc.m.|x) = V (101) 

L—¥(X) 

qGL* 

is a continuous nonnegative L-periodic function where 
^ H“ ‘ ‘ ‘ H“ — ^'i-j y 

and 

= E = i™,= 1- (102) 

qCL* 

We think this to be typical for crystals; see also the re¬ 
sult of the paper [5] for the ground state in ID and the 
forthcoming discussion in Section VH. 

(ii) Suppose that the extension of nk to (R*^)*, 

flL = E UkSw (103) 

keA* 

tends to a normalized regular distribution n; thus, 
J n{k) dk = 1. Then 

/i(x)= lim p-i(0|pi|x)= / h(k)e-*’" ’^ dk, (104) 

L—^cx) J 

a continuous nonnegative function decaying at infinity, 
and 

= I h(k)dk = ^im y(x) = 1. (105) 

This is typical for the single-particle momentum distribu¬ 
tion of bosons in the absence of Bose-Einstein condensa¬ 
tion. If there is BEC, in the limit we have h = hoi5o + fic 
with ho > 0 and he continuous, so that 

/^(x) ho as X —>■ oo. 

The non-decay of is the sign of the off-diagonal long- 
range order [21]. Periodic order can be seen also on the 
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limit of n^: if this contains a point measure at some 
k 0, in infinite space there is periodic order parallel 
to k [30, 31]. In all dimensions = 1 always holds, 
otherwise the kinetic energy density would diverge in the 
thermodynamic limit [31] . 

(iii) Suppose that tends to an unnormalized regular 
distribution v\ thus, 

0 < / f;(k)dk < 1. 


inequality estimates the distribution of the weight under 
the central peak. 

Lemma V.l. For any n > 2, 


1 

27r 



Fn (x) dx > 


- 4 -, if 712 < a < 27r 

V'J'TT 

if 0<a<7l2. 

(114) 


Proof. The result is obtained by integrating the lower 
bound 


Then for any x ^ 0, 

/-“•(x) < hm /-“ (y) = < 1, (106) 

y-+0 

cf. Eqs. (87) and (98), while /‘^ “ (O) = 1. Thus, /'=•“■ 
is discontinuous at zero, and the value of its jump is 
In Section VII we shall argue that the situation 
described in this point is characteristic to fluid phases. 
Outside the origin is a continuous function which 

in ordinary fluids decays at infinity. A special case is 

i> = i>oSo + Vc (107) 


with vq > 0 and Vc continuous, when 


^c.m.(x) 


f^(k)e-*’" ’^ dk ^ as X ^ oo. (108) 


We shall associate this case with superfluids. 


We now prepare some further results about r<oo = 1 
and Foo = 1- Combining Eqs. (80) and (81), 

Tl (^) = L-'^^/l(x) dx, (109) 

where 


Dn{x) 


l— — n 


sin[(2n + l)x/2] 
sin(x/2) 


( 110 ) 


Ftl{x) > max < n ( 1 — 


2 2 
n X 


12 


,07 □ 


(115) 


Next, we estimate the number of peaks that carry a 
prescribed weight. 

Lemma V.2. For any e > 0 there exists an integer Iq 
such that for any n > 21q, 


1 

27r 


^27rZo/^ 


Fn{x) dx > 1 — e. 


(116) 


> —27rZo / n 

Proof. Let I be any integer in the interval [l,n/2]. Then 

2»r(i + l) 


Fn(x) dx < 


1 


< 


1 


n2 [i_cos2^] - 7r[l-f|]Z2’ 


1 

27r 

''n L nJ L“ _ 

(117) 

Here we used 1 — cosx > x^/2 — x^/24 and l/n < 1/2. 
Because X)7i convergent, we can choose Iq such 

that 


1 


TT [l — 


12 J 


^r2<£/2. 


l=lo 


Thus, for any n > 2lo, 




(118) 


(119) 


is the Dirichlet kernel. A well-known fact is that 

f li:>„(x)l dx = O(logn). (Ill) 

J — TT 

We will also need the Fejer kernel, 

sin^(nx/2) 1 — cos(nx) 


Fn(x) = -Y = 


1=0 


nsi7(x/2) n(l-cosx)' 


( 112 ) 


Fn(x) is nonnegative and 

^ J F„(x)dx = l. (113) 


implying the result. □ 

Lemma V.3. Let be a sequence of real equicontin- 
uous and uniformly bounded functions in a domain 27 
of Define /"^(x) = limsupjgjr/(x) and /“(x) = 
liminf/g^/(x). Then are continuous in 27. 

Proof. We give the proof for /+. Fix an e > 0 and an 
X G 27. Because of equicontinuity, there exists a <5 > 0 
such that for any f £ F, |/(x) — /(y)] < e if Jx — yj < S. 
Let y be such a point of 27. Let {gn} and {h„} be two 
subsequences from F such that (/n(x) —>■ /“''(x), h„(y) —>■ 
/'’’(y) as n —>■ oo. There exists some Aq such that for 
any n > Ng, 

| 5 „(x) - /+(x)l < £, \hn{y) - /+(y)l < £, 


The graph of F^ consists of a sequence of peaks sepa¬ 
rated by zeros at 27r//n, I = ±1, ±2,_ The following 5n(y) < f'^iv) + ^n(x) < /''’(x) -|- £, 
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the last two inequalities following from the definition of 
the limsup. On the other hand, 

| 5 „(x) - 5 „(y)| < e, |/i„(x) - /i„(y)| < £. 

Thus, 

\9n{y) - /■^(x)l < 2e, |/i„(x) - /+(y)| < 2e. 

We may suppose that /'''(x) < Then 

/^(y) -2e < /i„(x) < /+(x) + £, 

that is, 

0 < /+(y) -/+(x) < 3£. □ 

We shall apply this lemma with = {/l}. Now 
/■(x) = liminf/L(x), /+(x) = limsup /l(x). (120) 

Because 0 < /l(x) < /l(0) = 1, we have 

0<r(x)</+(x)</±(0) = l. (121) 


only a part of those in the cube, the sums over the rect¬ 
angles and, hence, also their average is smaller than the 
sum over the cube. Substituting (/?L(k) from Eq. (80) and 
summing inside the integral, after a change of variables 
we obtain 



(124) 


where ||y|| = maxi|yi|. Due to equicontinuity and 
/l( 0) = 1, if £ > 0 is given then for any small enough 
A and any x such that ||x|| < A, /l(x) > 1 — £. Con¬ 
sider only L so large (or A so small) that L/X > 2. Let 
77 = 2nX/L, then 



s) 


1 

27r , 


1 d 


Fn{x) dx 


-V 


(125) 


Choose an integer Iq such that the inequality (118) holds 
true, and let n = n{L,e) = [loL/X\ > 2lo. Then 




/ 27m{L, e] 

1 ^ 


^ > (1 - £)‘^+b 


Theorem V.l. Suppose that the family of functions 
{/l} is equicontinuous in a neighborhood TX of the origin. 
Then /^(x) is continuous in V, so limx-i.o (x) = 1. 
Moreover, r<oo = 1, and for any sequence kl going to 
infinity, Tl{kl) 1- 

Remarks. 1. Instead of the equicontinuity of the whole 
sequence /l, equicontinuity of some subsequence 
would suffice. 2. By the Arzela-Ascoli theorem, the se¬ 
quence {/l} has a [not necessarily unique] continuous 
limit function in TX. 

Proof. The continuity of follows from Lemma 3. 
Equicontinuity implies that 

/-(x) < /(x) < /+(x) 

where we can use limsup or liminf in the definition of /. 
This, together with Eq. (92) yields r<oo = 1- Below we 
give also another proof for the saturation of E, that we 
shall use for the proof of ri(Ki) —>■ 1. Define 


Now 


therefore 


2TTn{L,e) 2 ttIq 
L - “ 



> lim gl 



>(i 


sY+\ 


so that 

limrf^^)=l. (126) 

£^0 \ X J 


Because E is monotone increasing and bounded by 1 (or 
because ?o/A tends to infinity as £ goes to zero), Eq. (126) 
is equivalent to limK->.oo E(k) = 1. Einally, consider any 
sequence kl going to infinity and choose n = \_Lhl/2'k\ . 
Then n > 21 q for L large enough, and 

lim El(kl) > (l-£)^+i 
L—¥00 


<JL 



then 


n—1 n—1 

( 122 ) 


El 


/ 27rn 

1“ 


> ctl 


(¥) 


(123) 


for any £ > 0; thus, the left member must be 1. □ 

In what follows, we discuss the relation between Eqo = 
1 and the properties of the sequence /l • 

Proposition V.4. If (/3 l(0) = o[L then 

limL-> 00 f — 0 almost everywhere, and Eoo — !■ 
Proof, (i) 


Indeed, TL{27rn/L) is obtained by summing in 

the central cube of A* of side length 27rn/L, while for 
CTL(27rn/L) we sum in centered rectangles inside this 
cube, and then average over the rectangles. Since all 
the terms are nonnegative and each rectangle contains 


(^l( 0) = L /L(x)dx = o(L 

JA 

implies that /^/L(x)dx = o(l). Then limL->oo /l (x) 
must vanish apart from a set of zero Lebesgue measure 
[which is nonempty since /l(0) = 1 for each L]. 
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(ii) The number of lattice points in A* such that ||k|| < k 
is not larger than (kL/tt) . Because for 

any k < co 

^l{k) < (kT/tt)'^ —>■ 0 (A —>■ oo) 

and, thus, r<oo = 0. □ 

Theorem V.2. Suppose that there exists some Kq > 0 
such that r(Ko) = 0- Then limL->.oo/ l(x) = 0 almost 
everywhere in 

Proof. By the monotonic increase of F, T{k) = 0 for 
every k < kq. Due to ri(K) > crL(K), c.f. Eq. (123), 

limsupcrL)^) = 0 if K < kq. 

L—¥00 

For notational simplicity assume that kL/2tt is an inte¬ 
ger. Because the integrand in Eq. (124) is nonnegative. 


By a variable transformation, Eq. (Ill) is equivalent to 

|^L|2 


~L/2 


D 


kL /27r 


/ 2'kx 
\~ 


dx = OilogKL), (131) 


so there is some constant C such that 


TUi^) <C{log kLY 



D^(^) 

dx 

kuU 

D^{^) 

dx 


(132) 

If L goes to infinity, the right member goes to zero, so 
lim^-j-oo Fl(k) = 0 for every k, which was the claim. □ 


VI. FREE ENERGY DENSITY 

Because the upper bound in Eq. (13) is subexponential 
in N, the free energy density can be obtained from Zi^red, 
c.f. Eq. (11). This we can rewrite as 


arin) > {2 tt) [ /l 

J\\y\\<27r^/KL 



d 

Y[P'<L/2^{yi)dy. 

2=1 


(127) 

The boundary value 2Y^ f kL of the domain of integration 
corresponds to a = tt in Eq. (114). It is easily seen that 

Fn{TTln) > dnliY, 


and Fn{x) increases monotonically when \x\ decreases 
from 7r/n to zero, where its value is n. Thus, in the 
domain of integration F)tL/ 27 r(yi) varies within a factor 
3. Applying Eq. (114) with a = tt, we obtain that 


(27r 


,—d 


T[F^L/2^{yi)dy > 

iiyii<2^v«i r=i 


In Eqs. (127) and (128) we substitute yi 
arrive at 



(128) 
(flixIV)Xi and 


(Jl{F) > 



/||x||<7r/Kni=l ■P'fi (^) dx 


(129) 

Fixing K < kq and letting L go to infinity, aL{K) goes to 
zero, so /l(x) must go to zero for almost every x in the 
cube {||x|| < tt/k}. Since this holds true for arbitrarily 
small K, the result follows. □ 

Note that the condition of the above theorem is fulfilled 
if Foo = 1- Therefore, we have the following. 


Corollary. In one dimension /‘^“ (O) = 1 and 

/c.m.(x) _ Q fQj. almost all (probably for all) x ^ 0 . 
Next, we prove a partial converse of Theorem V.2. 

Theorem V.3. Suppose that lim(log LY/^{x) = 0 
for almost every x G Then Fqo = 1. 

Proof. From Eq. (109), for n = k.L/2tt, 


rL(«) < jjL{^)\{ D.L/2. 


dx. (130) 


birred = ^ e-'^^<’<(q), (133) 

where 

fq = A-q.D ~ Fq q (134) 

and 

OO 

D^(q) = ^ (135) 

n=0 

Actually, it suffices to keep the q = 0 term of Zirred- Let 
qmax be the maximizer of Ai^(q) = —/3eq -|- lnr2^(q). 
Because contains N‘^ terms. 

The free energy density / = —limi_>oo L~ 
be given by 

/ = eo - lim T“‘'A:^(qmax) 

L—^oo 

where eo = limL“‘^Eo,o- From 

i/q = 

and Eq. (39) it follows that qmax = 0 in the case of 
Boltzmann or Bose statistics [32]. We conclude that 

f = eo-p-^ lim L-‘^lnD^(0). (139) 

L—¥00 

The free energy density is the Legendre transform of 
the microcanonical entropy density or, equivalently, 

/ = eo-k pmin {e(s) -/3“^s} , (140) 

where e(s) is the excitation energy (energy measured 
from Eq q) per particle at a (dimensionless) entropy per 
particle s in the thermodynamic limit. 


(136) 
In Z will 

(137) 

(138) 
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Proposition VI. 1. 


e(s) = e(0,s), (141) 


where e(0, s) is the limit of the excitation energy per par¬ 
ticle in the ensemble restricted to zero total momentum. 

Proof. We write 




pOO 

/ e-^^d 


Jo 

- 




(142) 

Here = 1 -|- max{n : E > — Ec^^o}, the 

microcanonical partition function in an ensemble where 
q is also fixed. The entropy SE^{E) increases with E, 
and its inverse E^^ (S) increases with S. Because the 
sum (135) converges for all /3 > 0, E^^ (S) must increase 
faster than linearly. Actually, if for a fixed q, S/N —>■ s 
as V —>• oo, then EE^{S)/N —>• e(q, s) which is a convex 
increasing function of s [33]. Therefore, the maximum of 
—/3e(q, s) -I- s is attained and 


_ gJV maXs>o{-/3e(q,s)-|-s}-|-o(JV) 


(143) 


Combining Eqs. (137)-(140) with Eq. (143) yields 
Eq. (141). □ 

The above result implies that on the level of thermody¬ 
namic functions the canonical ensemble is equivalent to 
the microcanonical ensemble restricted to Q = 0. Strong 
equivalence, meaning that the average value of local ob¬ 
servables is the same in the two ensembles, holds very 
probably as well. 

VII. TOTAL MOMENTUM AND THERMODY¬ 
NAMIC PHASES 

The examples given in Section V were meant to repre¬ 
sent the distribution of the total momentum in different 
types of thermodynamic phases. Our point was that e 
and 


/'=“ (x) = lim p l(0|pc.m.|x) 

L—¥C0 

are qualitatively different in fluid, solid, and superfluid 
phases. Below we formulate some detailed conjectures in 
this regard. The first two conjectures complete the theo¬ 
rems of Section II. Those results are valid both for bosons 
and fermions, with or without interaction. In noninter¬ 
acting systems something stronger can be expected. 

Conjecture VH.l. In two dimensional noninteracting 
systems for any /3 and p the asymptotic distribution of 
P/-\/]V is normal. 

In 3D systems of bosons we have to count with BEC. 
In the noninteracting case the fluctuations of Nq in the 
regime of condensation are of order [38]. Therefore, 
the following conjecture is slightly more challenging. 

Conjecture VH.2. In 3D, for noninteracting bosons, 
at all tempera tures and densities the asymptotic distri¬ 
bution of P/ \/N — (Nq) is normal. 


The verification or falsification of these two conjectures 
is within the reach of the existing mathematical meth¬ 
ods. The intuition behind them is that in the absence 
of interactions only an uncorrelated random motion of 
particles is possible and this must lead to a central limit 
theorem [39]. 

In what follows we focus on the interacting case. All 
the conjectures below refer to quantum systems with 
strongly tempered interactions. Strong temperedness is 
a condition on the decay rate of the interaction, nec¬ 
essary in order that the potential energy can be made 
periodic and that we can use periodic boundary condi¬ 
tions [33, 40]. 

Conjecture VH.3. In interacting systems above one 
dimension, at any temperature and density > 0. 

That is, in the infinite system the total momentum is 
finite with a nonvanishing probability in all the thermo¬ 
dynamic phases. 

The existence of collective excitations that carry a fi¬ 
nite momentum provides a strong experimental evidence 
supporting this conjecture. Such collective modes exist 
in all the thermodynamic phases; think about the den¬ 
sity waves in gases and liquids, the lattice vibrations in 
crystals. Strictly speaking, these excitations do not man¬ 
ifest themselves in thermal equilibrium. To see them, the 
equilibrium - specifically, the evenness [± symmetry] of 
the momentum distribution - must be broken, which is 
done by the experiment. ’Collective’ means that a macro¬ 
scopic number of particles contributes to the formation 
of each of these modes; thus, their wave vector is pro¬ 
portional to the total momentum which, therefore, must 
be finite with a positive probability. One dimension is 
an exception. Because in ID the probability that the to¬ 
tal momentum is finite tends to zero as the size of the 
system increases [T^™- = 0], we arrive at the following 
conclusion. 

Proposition VH.l. In one dimensional quantum sys¬ 
tems with strongly tempered interactions, at any positive 
temperature the density waves disappear in the thermo¬ 
dynamic limit. 

The restriction to T > 0 is crucial. The ground state 
can be periodically ordered [5], and there certainly exist 
low-lying excited states which correspond to its periodic 
modulations, i.e., to density waves. The absence of den¬ 
sity waves at T > 0 is the synonym of the absence of 
structural phase transitions. 

In the subsequent discussion of thermodynamic phases 
the use of states in infinite volume will be helpful. These 
are positive normalized linear functionals over the alge¬ 
bra of quasi-local observables. Here we introduce them 
through a strict minimum of complication. We suppose 
that all the limits below exist. Let 

OO OO 

zb = E i: 

||Q||<Kn=0 ||Q||>Kn=0 
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s<. = [zU E 

||Q||<Kn=0 


= [zi.Y E Ee”''''‘^’"nQ,n. 

||Q||>k n=0 
Then Z = Zh^ + and 


= 


^<K ^ ^>K 

which tends to 


^<K T ^>K 

yL ryh 

__ <^L I _ ^'>K _ qL 

rylj I rylj <K . 7^/^ '>K'i 

^<K "T ^>K 


(144) 


(145) 


5 = r---(s;) 5<« + [1 - r-“-('^)] ^>« 


(146) 


where 


Conjecture VH.4.: Gases and liquids. Above one 
dimension, in the fluid phases both probabilities, that of 
having a finite or an infinite total momentum in infinite 
space, are nonzero, i.e., 0 < < 1. Density waves 

have a finite total momentum whose probability distri¬ 
bution is continuous. The infinite total momentum is the 
resultant of random single or few-particle motions which 
in finite volumes create a total momentum of the order of 
\/N. Thus, in Eq. (152), i> is continuous and > 0. 

The difference between gas and liquid is only quantita¬ 
tive. When cooling the system, it is possible to go around 
the critical point by continuously varying the density. 
When cooling at a constant [global] density, there will 
be a phase transition with phase separation: the liquid 
phase has a higher density and a higher probability of a 
finite total momentum than the coexisting gas phase. 


S = lim 5^, 5 <k = lim 5<^, S'>k = lim S'>^. 

L —¥00 L —¥00 ~ L —¥00 

(147) 

These limits are taken on local observables. The decom¬ 
position can be done for any k. Sending k to infinity, 
with 


5<oo = lim S<K, Soo = lim 5 >k (148) 

K—^OO K—^OO 

we have 


5 = r-“-5<oo + r-“-5^. (149) 

The first term can further be decomposed. It reads 

r<“- 5 <oo = J dk, (150) 


where 


(j(k) = lim lim 
g—>-0 L—^00 


El 


Q-k||<g 


y^oo „-I3Eq,„yt 

Z^n=0 ^ 


El 


Q— Z^n—0 


Z^n =0 ^ 


(151) 


The final form of S is 

5 = y’:>(k)a(k)dk + r^“,Soo. (152) 

In Eq. (152) the first term is composed of finite momen¬ 
tum states and represents the collective modes; the sec¬ 
ond term is composed of infinite momentum states and 
represents the few-particle excitations. 

Few-particle motion is a coordinated motion of a small 
group of particles which evolves independently from the 
rest of the system. Rotons are thought to belong to this 
category [41]. In 3D at high energies, in 2D at all energies 
such nearly separated small groups occur with a positive 
density in a non-vanishing fraction of eigenstates. The 
sum of the group momenta will typically be of the or¬ 
der of ■\/N. There exist also excited states whose total 
momentum is finite and is the sum of a finite number of 
single- or few-particle momenta. However, at T > 0 the 
statistical weight of these states becomes negligible with 
the increasing system size. 


Conjecture VII.5.: Solids. In an infinite solid phase 
the random few-particle motion vanishes almost surely, 
implying that the total momentum is finite with proba¬ 
bility 1. In a crystal the probability distribution of the 
total momentum in infinite space is concentrated on a 
lattice, cf. Eq. (100). So = 0, f> = X^keL* 

/ iz(k) dk = 1 , and 


S = E] ^'k'7(k). (153) 

keL* 

In view of what we found in 2D, this conjecture ex¬ 
cludes the solid phases in 2D at any positive tempera¬ 
ture. Indeed, by Theorem II.2 and Proposition V.3 we 
know that T^”' > 0- This, together with Conjectures 
VII.3 and 4 predicts that 2D systems exist only in fluid 
[gas, liquid or superfluid] phases. This conclusion is in ac¬ 
cordance with the preservation of translation invariance 
in 2D particle systems in continuous space [42-45]. In 
addition, there is the interesting consequence that if the 
hexatic phase were to exist classically, quantum fluctua¬ 
tions would destroy it: Orientational long-range or even 
quasi-long-range order is incompatible with an isotropic 
random few-particle motion which implies exponentially 
decaying correlations. 

A special kind of crystals are the so-called coherent 
quantum crystals, introduced in the papers [46, 47]; see 
also Ref. [31]. They are thought to be somewhere be¬ 
tween crystals and fluids, showing periodic long-range 
order, EEC and fluid properties (large kinetic energy) si¬ 
multaneously. In our interpretation their simplest repre¬ 
sentative would be characterized by an asymptotic prob¬ 
ability measure of the form 

+ i>q[Sq + (5-q] + J f^c(k) dk -F (154) 


where q 7 ^ 0, f>o > izq > 0, and bo -|- 2bq < 1. Bosons in¬ 
teracting with an integrable pair potential whose Fourier 
transform has a negative minimum can be in such a state 
at very low temperatures; q would be close to the mini¬ 
mum. 
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Conjecture VII. 6.: Superfluids. A superfluid is a 
fluid in which the asymptotic probability that Q = 0 is 
nonvanishing. That is, lim^^-j-o r'^ ™ (K) = > 0. 

Note that limi_>oo may vanish while t'o > 0. Ac¬ 
cording to this definition, a positive t'o should appear in 
2D systems of bosons [helium films, trapped Bose gases] 
and in 3D liquid "^He when the temperature decreases 
below the transition temperature Tg. We remark that 
do > 0 also in crystals [{dkjkGL* is positive definite, 
therefore vq > dkj; the specificity of superfluids com¬ 
pared to crystals is that outside k = 0 the asymptotic 
distribution v is continuous and > 0. This implies 
that has a jump at 0 and tends to the constant 

do > 0 as X —>■ oo, without exhibiting any long-range 
order. So we have 

do+y’d,(k)dk-|-r-“- = 1, (155) 

all the three terms are positive, and dc is a continuous 
measure, cf. Eq. (107). The two-fluid state is 

5 = do 5s -I- [1 - do] 5„, (156) 

the super and normal fluid states are 


Then H'' and P'' are, respectively, the energy and the 
total momentum operators in the comoving frame [48]. 
H'' has the same eigenstates as H. Because mvlh G A*, 
the eigenvalues of iF'" are also unchanged [so Tr = 

Z], only their assignment is permuted: 

^ ‘^Cl+Nmv/hyn ~ ‘^Cl+Nmv/h,n: 

^ ‘^Cl+Nmv/h^n ~ ‘^Cl+Nmv/h^n- (160) 

Moreover, 

(161) 

Equations (160) and (161) together with Eq. (8) express 
Galilean invariance [5]. The density matrix reduced to 
the center of mass becomes 

Pc.m. = ^ X! ’^q\Q + Nrnv/h){Q + Nmv/h\, (162) 
qga* 

and therefore 

(0|Pc.m.|x) = e-^'''’'(0|pc.m.|x). (163) 


5s 


Sn 


a{0) 

[1-do]-^ 


J dc(k)CT(k) dk-|-r[] 


‘•5. 


(157) 


Thus, the super and normal fluid components are not 
parts of the system but its alternative states that com¬ 
pose the infinite volume Gibbs state by a convex com¬ 
bination. Gonjecture VII.2 implies that do = 0 in the 
noninteracting Bose gas which is, therefore, not a super¬ 
fluid [39]. 

This image of a superfluid is radically different in many 
respects from what the well-established theory offers to 
us. One point is particularly important: in Eq. (156) 
the superfluid state has the same temperature and the 
same energy and entropy density [cf. Proposition VI. 1] 
as the normal fluid state, and not zero temperature and 
zero entropy density, as Landau suggested it [6]. This 
is in accordance with the fact that the A-transition is 
an equilibrium phase transition. In the discussion of the 
dissipative flow below we shall return to this question. 

The above equations describe the superfluid at rest. 
Now we consider the superfluid in motion. A symmetric 
V-particle wave function $ describes a system in motion 
of velocity v with respect to its container A if 

($|pj|$)=mv (anyj). (158) 

Because then (d’jpj — mvj$) = 0, in the comoving refer¬ 
ence frame [the one moving together with the system] the 
single-particle momenta will be pj = Pj — mv. Choose 
V so that mv(h G A*. Since A* fills densely the whole 
space when L —>■ oo, this is not an important restriction. 
Introduce 

1 ^ 

= + = (159) 

i=i i 


Its square root can play the role of a macroscopic wave 
function, 

^(x) = \/(0lpXm.|x) = e"^'''^\/(0lpc.m.|x). (164) 

Eor large N and large x, 4'(x) ~ Because 

the function under the square root is nonnegative, by 
taking the positive square root the phase of is 

N 7Tl 

Arg4'(x) = —^v-x, (165) 

so that the velocity is obtained as 

V =-^^Arg4'(x). (166) 

The construction (144)-(152) of the states in infi¬ 
nite volume can be repeated by replacing IIq „ with 
^Q+Nmv/h,n but leaving unchanged. All the 

infinite-volume states will be labeled by an upper index 
V. The infinite volume Gibbs state (152) becomes 

5- = I ^(k)a''(k) dk + r---5-. (167) 

The two-fluid state of the superfluid is 

5" = :>o5J + [l-f>o]5)[, (168) 


with 


5: 

5: 


= (^"(0) 

= [1 - i>o] 


-1 


j>,(k)a-(k)dk-kr---5^ 


(169) 
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and unchanged super and normal fluid fractions. In gen¬ 
eral, V remains unchanged because it depends only on the 
energies. For the same reason, the free energy density of 
S'^ agrees with that of S. 

The situation described above corresponds to a hypo¬ 
thetic frictionless flow. In a realistic flow the velocity de¬ 
pendence of the density matrix is not Galilean invariant. 
Below we show how dissipation leads to a critical velocity. 
Imagine that a constant flow velocity v is maintained e.g. 
by moving the capillary with velocity —v, and v = |v| is 
smaller than the critical velocity Ucr(T) that we are going 
to determine. Now v will be the velocity only in the su¬ 
perfluid state, Vs = V. In the normal fluid state friction 
on the walls and viscosity reduce the mean drift veloc¬ 
ity to v„, |v„| < V. The accompanying loss of kinetic 
energy is converted into heat which partly is dissipated 
into the environment, partly raises the temperature in 
both the super and the normal fluid states. In effect, 
the temperature of the superfluid component cannot be 
lower than that of the normal component, otherwise the 
superfluid would immediately disappear! This point is 
important because it contradicts Landau’s intuition [ 6 ], 
therefore we spell it out. Suppose that in the superfluid 
state the temperature is T, in the normal fluid state it is 
T', and T <T' <Ts. Then in finite volumes the density 
matrix (145) becomes 


where 77 < 1 is the efficiency, i.e., the fraction of the heat 
that raises the temperature. Now 77 may also depend 
on t, but we suppose that it does not tend to zero, and 
neglect this dependence. In this case a{t) —>■ 0 as t —>■ 00 . 
At time t the state of the system is 

SLne^iTt) = i>o{Tt) SJ{Tt) + [1 - Pom sm^'m- 

(174) 

The state is out of equilibrium, because the two terms are 
generated by different Hamiltonians, and j.g_ 

spectively. The temperature Tt increases in time; if v was 
small enough then limt_>oo Tt = Ty < Tg, and 5^oneq(^i) 
converges to a non-Gibbsian steady state whose normal 
component is at rest in A while the superfluid flows with 
velocity v: 

SZone^m = 5J(r„)+[l - Po{Ty)] 5„(T„). (175) 

From Eq. (173), 

Ty = T + r]mv^/ks- (176) 

As V increases, Ty tends to Tg, Po{Ty) goes to zero, and 
thermal equilibrium is restored: the system will be in 
the state 5„(Ts), at rest with respect to A. The critical 
velocity is obtained from the equation = Tg which 
yields 




+ 


zUt) 

{t) + zUt') 

z^m)+ zut') 


■ skm) 
sUt')- 


(170) 


The superfluid/normal fluid ratio is the double limit of 

ZLiT) ^ ZUT) ZUT) .^ 7 ^. 

Z^T') ZUT) ZUT')' 

Applying lim^-^o limi_>oo to this equation, the first frac¬ 
tion on the right remains positive because T < Tg. How¬ 
ever, the second fraction is exponentially small in N and, 
thus, goes to zero as L goes to infinity. [In finite volumes 
the ratio of two partition functions of the same system 
taken at different temperatures is exponential in N.] 

The velocity of the normal fluid is time dependent, 
Vn = v„(t). Supposing that v„(0) = v, we can write it 
in the form 


v„ = a{t)v, a(0) = 1. (172) 

As long as v„ ^ 0, friction is in action, therefore v„ de¬ 
creases and the heat production continues. Let Tt denote 
the temperature of the system at time t, Tq = T. Since 
only one degree of freedom, the single-particle momen¬ 
tum parallel to v is concerned, conservation of energy 
dictates 


^ksTt - ^ksT = - aitf], (173) 


Vcr{T) = y^kB{Tg-T)/{7jm). (177) 

A possible definition of the ground state in infinite vol¬ 
ume is limT-s-o ^(r). A critical velocity for this state is 
provided by Ucr(O) = ^/ksTg/{'qm). For He H at satu¬ 
rated vapor pressure Tg = 2.17 K, which gives Ucr(O) = 
(67/v^) m/s > 67 m/s. Because of the walls, in a capil¬ 
lary there is normal fluid and dissipation also at T = 0, 
so the above discussion is relevant: the dissipation will 
heat the system to the temperature Ty = qmv^ /ks- 
The mechanism by which friction affects the moving 
normal fluid but not the superfluid can be understood 
qualitatively. In the moving system friction on the walls 
is a surface effect, internal friction (viscosity) is a bulk 
effect, their best description is via a spatially random 
uncorrelated perturbation. Such a perturbation will not 
excite collective modes, it can influence only small sep¬ 
arated groups of particles and excite them to a (macro- 
scopically) higher energy. These groups then relax and 
emit incoherent radiation which partly leaves the sys¬ 
tem, partly is reabsorbed and heats thereby. Now 
must be an infinitely entangled state, with no separable 
small groups of particles in it. Separated small groups 
of particles can be found in a macroscopic number only 
in the eigenstates 'tpQ+Nmv/h.n = with 

IQI ~ VN, which give rise to the component of S^- 
Our definition of the superfluid fraction Pq is in obvious 
analogy with the condensate fraction ho. The connection 
between superfluidity and BEG is a subtle question. In 
two dimensions there is no BEG [16-18] but superfluid¬ 
ity does exist [19, 20]. On the other hand, even though 
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only an increased peak height and not a sharp peak could 
be measured at k = 0 in the single-particle momentum 
distribution [24], there is today a consensus based on an¬ 
alytic arguments, experiment and numerics [21-26] that 
in bulk liquid helium the superfluid transition and BEC 
occur simultaneously. However, the condensate ratio sat¬ 
urates at less than 10% while the whole liquid becomes 
superfluid as the temperature goes to zero. So one cannot 
simply identify the superfluid with the Bose-condensate: 
apart from the numerical mismatch, they are conceptu¬ 
ally different, the condensate being a part of the sys¬ 
tem while the superfluid is, according to us, a state of 
it. According to the Bogoliubov-Landau theory, in three 
dimensions BEC acts as a catalyst that triggers the su¬ 
perfluid transition. However, it may well be the other 
way around: when two phenomena are simultaneous, it 
is difficult to decide, which one is the cause and which is 
the consequence. We prefer to consider the emergence of 
the Q = 0 subspace as the primary event that causes the 
superfluid transition both in 2D and in 3D and, addition¬ 
ally, BEC in three dimensions. Among many other signs, 
this also shows that the mechanism of BEC for interact¬ 
ing bosons must be quite different from the saturation 
effect which is at the origin of the condensation in the 
three dimensional noninteracting Bose gas. 

VIII. COMMENT ON LANDAU’S CRITERION 
OE SUPERELUIDITY 

The physical source of the critical velocity Vcr{T), de¬ 
rived in the preceding section for superfluids in a capil¬ 
lary, is quite different from that of the critical velocity in 
Landau’s theory. By analyzing Landau’s original publi¬ 
cation [6], we show that his argument about the condition 
of excitability of the superfluid ground state in a dissipa¬ 
tive flow contains an error, which consists in deriving the 
critical velocity from a frictionless flow. Indeed, he ob¬ 
tained the critical velocity through a mere Galilean boost 
which, by preserving the energy spectrum, does not lead 
to dissipation. 

The details are as follows. In Section 4 of Ref. [6], Lan¬ 
dau investigated the stability of the superfiuid ground 
state of helium liquid in a capillary against low-energy 
excitations that a flow of velocity v can create by losing 
kinetic energy. His argument is based on the properties 
of eq = Eq^o — Eq^o, the energy gap to the lowest-lying 
state of momentum HQ, that we introduced in Eq. (134). 
This quantity often appears in rigorous works about Bose 
systems, see e.g. Refs. [34, 35]. In He H, eg is measured 
up to about 4A“^. Its qualitative features, known al¬ 
ready to Landau, very probably do not depend on T, 
provided that T < T^: the curve starts linearly, passes 
over a maximum and exhibits the famous roton mini¬ 
mum [10, 11, 36]. If eg could be measured at T = 0, it 
would certainly show the same features as in the mea¬ 
surement [36], done at 1.1 K. In our notations, Landau 
wrote down the equation 

eq+jVmv/?i = Cq -k ftq • V -h ^Nmv'^ , (178) 


first with Cq « c/i|q| near q = 0 [Eq. (4,1)] and second, 
with Cq « A -|- (2/r)“^fi^(q — q^.)^ [Eq. (4,3)] near the 
roton minimum q^.. He considered Eq. (178) as an equa¬ 
tion of energy balance for the moving fluid, and said that 
for a loss of kinetic energy and accompanying excitation, 

Cq -t- fiq • V < 0 or v > for some q (179) 

^|q| 

must hold. In particular, a phonon of momentum hq can 
be excited if 

c|q| -I- q • V < 0, (180) 

and a roton can be excited if 

A-I-(2^)“^li^(q — qr)^ — ?i|q|u < 0. (181) 

This, however, is not true, these inequalities do not point 
to any distinguished value of v. If we add Eq q to both 
sides of Eq. (178), we obtain 

Eq+Nm^/h,o = Eq^o + hc[ ■ V + (182) 

which is Eq. (8) with n = 0 and k = mv/h if we sup¬ 
pose, as everywhere in this paper, that mv/h G A*. By 
Galilean invariance [5], the set equalities 

{A'q.O -(- fiq • V -|- — VmV^lqgA* = {£'q-|-JVmv//l,o}qeA* 

= {A-q.olqSA*, (183) 

hold true: the first equality is term-by-term, the second 
is true with a permutation of the elements, which changes 
as V varies. In particular, as seen from Eq. (160), 

V’AT^v/fi.o = (184) 

the lowest-lying eigenstate that describes the system 
moving with velocity v with respect to A [i/^jvmv/n,o sat¬ 
isfies Eq. (158)], and V’o.O) the ground state at v = 
0 , exchange energy: that of the former decreases by 
Atov^/ 2 and becomes Eq^o, that of the latter increases 
by Atov^/ 2 and becomes Eq^o + Atov^/ 2. This occurs 
for any velocity such that mv/h G A*; the sound veloc¬ 
ity or the slope of the straight line drawn to the roton 
minimum are not distinguished. 

More generally, one observes that eg plays no role: the 
excitation energy supplied by the kinetic energy loss of 
the real flow [which is dissipative even at T = 0 because 
of the walls of the capillary] is of order N, highly above 
eg from |Q| = 0 to |Q| ~ '/N- Thus, if the critical 
velocity should come from a comparison of eg with the 
available excitation energy, then it would be zero. Note, 
however, that condition (179) is relevant when bulk He H 
is locally perturbed, by moving in it a small object with 
velocity v. 

eg is nonetheless good for one thing. At T > 0, still 
below the transition temperature, the zero point of the 
IQ] —>■ eg curve does not belong to the ground state. 
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but to a state of zero total momentum whose energy is 
in the order of N higher than that of the ground state. 
This state cannot be else than the superfluid state at the 
given temperature, providing an experimental support to 
our view, that the superfluid state is a state of zero total 
momentum. 

IX. SUMMARY 

In this paper we studied the probability distribution of 
the total momentum of quantum fluids at positive tem¬ 
peratures. Sections II-VI contain mathematical results 
about this quantity in finite volumes and in the limit 
of infinite space. The motivation behind this work is 
our conviction that the distribution of the total momen¬ 
tum is significantly different in different thermodynamic 
phases, and therefore can be used to characterize fluids, 
solids and superfluids. This conviction is expressed in 
the examples given in Section V and through the conjec¬ 
tures formulated in Section VII. The main assumption is 
that in interacting quantum systems above one dimension 
the total momentum remains finite with a nonvanishing 
probability in the thermodynamic limit; for a justifica¬ 
tion we referred to the existence of collective excitations. 
One may then envisage that with some probability the 


infinite system can be in a state of zero total momentum. 
This is precisely what we think to happen during crystal¬ 
lization and in the superfluid phase transition. The dif¬ 
ference between the two cases is that in crystals the total 
momentum is finite and is distributed on a lattice while 
in a superfluid it can be infinite and otherwise its distri¬ 
bution is continuous outside the origin. A measurement 
can project the system into any component of the ther¬ 
mal equilibrium state which has a non-vanishing proba¬ 
bility, such as the superfluid component or a state with 
a distribution over nonzero finite momenta representing 
phonons in normal fluids or in crystals. We suggested 
the use of the density matrix reduced to the center of 
mass for the definition of a macroscopic wave function, 
and commented on the effect of dissipation. Section VIII 
was devoted to a critical discussion of Landau’s criterion 
of superfluidity. The relation between BEC and the su¬ 
perfluid transition remains to be clarified. The major 
unsolved problem is, however, to prove the phase transi¬ 
tions. 
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